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Page  k,  line  5:  Replace  the  quantity  M^  -  1  by  \JM^  -  1. 

Page  11,  equation  (if)  In  the  first  value  for  p,  for  plan  form  D, 

T  7 

the  term  —  within  the  final  parentheses  should  be  — . 

2b  2b 

Page  21:  Correct  equations  (B3)  as  follows: 

In  line  2,  the  expression  should  be 

In  lines  4,  6 ,  7,  13,  15,  and  18,  all  primes  appearing  on  the  right- 
hand  side  of  these  equations  should  be  replaced  by  the  superscript  1. 
In  line  l4 ,  the  sixth  term  of  the  equation  for  -2rtM^ '  should  be 
5 3^12°  instead  of  5jHq2- 

Page  22:  In  line  5*  the  "Mq^"  within  the  bracketed  expression  should 
be  "M2.” 

In  line  7*  The  right-hand  side  of  the  equation  for  p^  should  be 
preceded  by  a  minus  sign. 

Page  23:  In  lines  2,  3j  and  4,  the  equations  for  and  5g 

should  be  corrected  as  follows : 
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Page  23j  next  to  last  line:  The  superscript  "vp"  should  be  the  super¬ 
script  "1/2." 

Page  25,  line  4:  The  equation  for  should  be  corrected  as  follows: 
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SUMMARY 


Expressions  for  the  velocity  potential  and  the  lift  and  moment  are 
derived  for  thin,  oscillating,  "arrowhead”  type  wings  with  supersonic 
edges  moving  at  a  constant  supersonic  speed.  The  triangular  wing  and 
any  one  of  the  series  of  wings  obtained  from  it  by  cutting  its  rearward 
part  so  that  the  flow  normal  to  the  resulting  trailing  edges  remains 
supersonic  are  included  in  the  term  arrowhead  type.  Explicit  results 
for  the  lift  and  moment  are  given  for  the  wings  undergoing  harmonic 
pitching  and  vertical  oscillations. 

Closed  expressions  for  the  velocity  potential,  section  force  and 
moment  coefficients  for  any  arrowhead  wing,  and  total  force  and  moment 
coefficients  for  only  the  triangular  or  delta  wing  are  developed  explic¬ 
itly  to  the  third  power  of  the  frequency  of  oscillation.  These  expres¬ 
sions  may  be  sufficient  for  most  practical  applications. 

The  wings  considered  are  found  to  exhibit  the  possibility  of  nega¬ 
tive  damping  in  torsion  for  certain  ranges  of  Mach  number  and  axis-of- 
rotation  position.  A  figure  showing  the  ranges  of  these  parameters  for 
triangular  wings  is  given  herein. 


INTRODUCTION 


In  high- speed-aircraft  design,  a  knowledge  of  the  air  forces  and 
moments  that  act  on  various  types  of  oscillating-wing  plan  forms  is 
often  required.  Such  knowledge  is  useful  in  considering  a  number  of 
instability  problems,  among  which  are  wing  flutter  and  low-frequency 
instability  of  the  aircraft  involving  control- surface  deflections.  The 
main  source  of  theoretical  information  has  been  the  solution  of  the 
linearized  differential  equation  for  compressible  flow. 

The  problem  of  finding  the  air  forces  on  an  oscillating  or  steady 
thin  wing  in  supersonic  flow  is  usually  formulated  in  terms  of  the  dis¬ 
turbance  velocity  potential  for  the  wing.  The  disturbance  velocity 


2 


MCA  TN  2k9k 


potential  is  constructed  by  superposing  solutions  of  the  linearized 
differential  equation,  corresponding  to  sources,  doublets,  or  higher- 
order  singularities,  in  such  a  way  as  to  satisfy  the  boundary  conditions 
at  the  wing  surface  and  at  infinity.  The  linearized  differential  equa¬ 
tion  and  its  boundary  conditions  constitute  the  boundary- value  problem 
for  the  velocity  potential  for  the  wing. 

The  present  paper  is  concerned  with  the  boundary- value  problem  for 
the  velocity  potential  for  an  apex-forward,  thin,  flat,  oscillating, 
triangular  or  related  wing  with  supersonic  edges  moving  at  a  constant 
supersonic  speed.  According  to  reference  1,  this  problem  may  be  clas¬ 
sified  as  "purely  supersonic,"  since  the  upper  and  lower  surfaces  of  the 
wing  can  be  regarded  as  acting  independently  of  one  another,  and  is  sat¬ 
isfied  by  a  surface  distribution  of  sources  with  local  strength  propor¬ 
tional  to  the  local  prescribed  normal  velocity  at  the  wing  surface. 
Although  reference  1  does  not  explicity  treat  the  oscillating  triangular 
wing,  its  general  solution  must  be  regarded  as  giving  the  velocity 
potential  for  this  wing  in  integral  form.  The  purpose  of  the  present 
paper  is  then  to  obtain  an  integrated  form  for  the  velocity  potential, 
from  the  integral-form  solution  of  reference  1,  for  the  particular  case 
of  the  oscillating  triangular  or  related  wing. 

In  a  strict  sense  this  integrated  form  applies  to  a  semi-infinite 
triangular  plan  form.  But,  since  the  wing  wake  is  of  no  concern  in  the 
present  development,  the  velocity  potential  may  be  considered  to  apply 
to  any  one  of  the  series  of  wings  obtained  from  a  triangular  wing  with 
supersonic  edges  by  cutting  the  rearward  part  of  this  wing  in  such  a 
manner  that  the  flow  normal  to  the  resulting  trailing  edges  remains 
supersonic.  All  these  wings,  including  the  triangular  wing,  are  hence¬ 
forth  referred  to  in  general  as  arrowhead  wings. 

There  are  several  other  papers  closely  associated  with  the  subject 
being  considered.  In  references  2  and  3  the  treatment  of  oscillating, 
finite,  swept  wings  involved  an  integral  form  for  the  velocity  potential 
quite  similar  to  that  for  the  oscillating  arrowhead  wing.  In  refer¬ 
ences  *4-  and  5,  more  directly,  expressions  for  the  total  lift  and  moment 
were  given  for  the  particular  oscillating  arrowhead  wing  with  super¬ 
sonic  edges  known  as  the  "wide"  delta  wing. 

The  present  paper  effectively  includes  the  total  forces  and  moments 
for  the  wide  delta  given  in  references  k  and  5-  In  addition  the'  present 
paper  gives  the  forces  and  moments  on  any  streamwise  wing  section  for 
not  only  the  wide  delta  but  also  for  the  more  general  arrowhead  wing. 

The  section  rather  than  the  total  forces  and  moments  are  desirable,  for 
example,  in  a  strip  flutter  analysis  that  includes  the  wing-flutter  mode 
shape. 

The  velocity  potential  for  the  arrowhead  wing  does  not  appear  to 
be  obtainable  in  terms  of  elementary  functions  (see  references  1  and  k) . 
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The  solution  of  this  type  of  problem,  therefore,  is  approximated  by 
various  kinds  of  expansion.  Thus  in  reference  4  the  velocity  potential 
was  expanded  in  a  series  of  Bessel  functions;  whereas  in  reference  6, 
in  the  treatment  of  the  rectangular  wing,  an  expansion  in  powers  of  the 
frequency  of  oscillation  was  employed.  The  frequency-expansion  method 
is  utilized  in  this  paper  and  only  the  first  few  terms  of  the  velocity 
potential  and  subsequent  forces  and  moments  are  obtained.  The  first 
few  terms  are  considered  adequate  for  a  large  class  of  practical  appli¬ 
cations.  In  the  derivation  of  the  velocity  potential  and  the  forces 
and  moments,  the  wings  are  assumed  to  be  undergoing  harmonic  oscilla¬ 
tions  in  pitch  and  vertical  translation. 


SYMBOLS 


w(x,y,t) 


disturbance-velocity  potential 
rectangular  Cartesian  coordinates 
velocity  of  main  stream 
velocity  of  sound  in  main  stream 
time 

normal  velocity  at  surface  of  wing  at  point  (x,y) 

function  defining  vertical  displacement  of  point  (x,y) 
of  wing 

abscissa  of  axis  of  rotation  of  wing  section  as  shown  in 
figure  1 

air  density  in  main  stream 

vertical  displacement  of  axis  of  rotation 

maximum  amplitude  of  vertical  displacement  of  axis  of 
rotation,  positive  downward 

angle  of  attack 

maximum  amplitude  of  angular  displacement  about  axis  of 
rotation,  positive  leading  edge  up 

time  derivatives  of  h  and  a,  respectively 
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cd  circular  frequency  of  oscillation 

M  free -stream  Mach  number  (u/a) 

p  =V'm2  -  1 

£,tj  rectangular  Cartesian  coordinates  used  to  represent 

location  of  sources  in  xy-plane 

up 

an,bm  functions  of  cd,  x,  and  M  defined  in  equation  (ll) 

X  slope  of  leading  edge  of  wing  as  shown  in  figure  1 


■a 


lim  (j) 

px — >1 


u,v  characteristic  coordinates  defined  in  appendix  A 


1  -  px 
7  “  1  +  px 

a  =  p2X2  -  1 


b  semichord  of  midspan  wing  section 

|i,V,I-Lq  nondimens ional  coordinates 

k  reduced  frequency  (bcn/u) 


Ap 


local- surface  pressure  difference 


P 

Ma 


Li, Mi 


M. 


a 


section  force  per  unit  span  on  wing  strip  parallel  to  main 
stream,  positive  downward 

section  moment  per  unit  span  on  wing  strip  parallel  to  main 
stream,  taken  about  axis  of  rotation  xq,  positive  leading 
edge  up 

components  of  section  force  and  moment  coefficients, 
respectively,  defined  in  equation  (Bl) ;  i  =  1,2,3,^ 

total  force  on  delta  wing,  positive  downward 

total  moment  on  delta  wing  about  axis  of  rotation  x  =  Xq, 
positive  leading  edge  up 


NACA  TN  2494 


components  of  total  force  and  moment  coefficients, 
respectively,  defined  in  equation  (B4);  i  =  1,2, 3, 4 


G  ,H  I 
mrr  mn 

l  functions 

of 

P,  and 

ai>pi^j 

7i»8iJ 

functions 

of 

M, 

and 

distance  to  point  at  which  extended  trailing  edge  of  plan 
form  D,  figure  3>  intersects  x-axis 


semi span  of  wing 


FORMULATION  OF  PROBLEM 


Consider  uniform  supersonic  flow  past  a  thin,  flat,  arrowhead  type 
wing  as  shown  in  figure  1  (with  its  leading  edges  outside  the  Mach  cone 
generated  by  the  nose  of  the  wing) .  The  wing  is  referred  to  a  coordinate 
system  fixed  in  space  and  is  assumed  to  he  creating  small  disturbances 
in  the  main  stream  flowing  past.  Then,  if  the  undisturbed  stream 
velocity  is  in  the  direction  of  the  positive  x-axis  and  is  of  magni¬ 
tude  U,  the  equation  satisfied  by  the  disturbance  velocity  potential 
for  the  wing  is 


±fd 

a2\3t 


(1) 


The  wing  is  assumed  to  be  performing  oscillations  of  small  ampli¬ 
tude  about  its  own  undisturbed  position,  the  plane  Z  =  0.  The  boundary 
condition  to  be  satisfied  at  the  wing  surface  is  then 


)g — ^0 


w(x,y,t)  =  U 


&Z„  dz. 


m 


m 


(2) 


where  7^  is  the  vertical  displacement  of  a  point  (x,y)  of  the  wing. 

Note  that  this  boundary  condition  is  evaluated,  in  accordance  with 
small-disturbance  linearized  theory,  at  the  xy-projection  of  the  wing 
rather  than  at  the  wing  itself.  The  additional  boundary  conditions, 
that  only  the  wing  can  support  a  pressure  difference  and  that  the  sources 
of  disturbance  must  not  be  felt  ahead  of  their  respective  Mach  cones, 
are  automatically  satisfied  by  the  type  of  source  synthesis  to  be 
employed  in  the  solutions. 
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After  the  boundary- value  problem  is  solved  for  the  velocity 
potential  (f),  the  pressure  on  either  wing  surface  may  be  found  by  means 
of  a  linearized  form  of  Bernoulli !s  equation 

+  u^)  (3) 

where  p  is  the  density  in  the  undisturbed  stream. 

For  the  present  problem  of  a  wing  performing  small  harmonic  tor¬ 
sional  oscillations  of  amplitude  <Xq  about  some  spanwise  axis  xq 
and  small  harmonic  vertical  translations  of  amplitude  hg  (see 
fig.  1(b)),  the  equation  for  is 

7^  =  h  +  (x  -  x0)a  = 

where  to  is  the  frequency  of  oscillation.  Substitution  of  equation  (k) 
into  equation  (2)  gives 

w(x,y,t)  =  h  +  Ua  +  (x  -  xQJa  (5) 

Since  equation  (l)  is  linear,  the  velocity  potential  satisfying  equa¬ 
tions  (l)  and  (5)  may  he  regarded  as  the  sum  of  three  potentials,  each 
being  associated  with  one  of  the  terms  on  the  right-hand  side  of  equa¬ 
tion  (5) •  Thus  the  potential  is  obtained  in  the  form 

*  -  K +  K  *  *k  <«> 


hp  *  (x  -  Xojap 


imt 


(M 


Thickness  effects  are  considered  negligible  and,  as  a  result, 
the  velocity  potentials  in  equation  (6)  are  associated  only  with  lift 
and  are  anti symmetrical  with  respect  to  the  plane  Z  =  0.  Thus  only 
one  surface  of  the  projected  wing  need  be  considered.  The  top  surface 
(Z  =  +0 )  is  chosen  in  this  analysis  and,  since  the  anti symmetrical 
potential  is  simply  opposite  in  sign  on  the  bottom  surface,  the  pres¬ 
sure  difference  supported  by  the  wing  is  obtained  by  means  of  equa¬ 
tion  (3)  as 


(T) 
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EQUATIONS  FOR  VELOCITY  POTENTIAL 
Integral  Form 


The  boundary-value  problem  for  the  arrowhead  wing  with  supersonic 
edges,  discussed  previously,  is  similar  to  that  for  the  wing  of  infinite 
span  treated  in  reference  1.  The  problem  in  reference  1  is  shown  to  be 
satisfied  by  a  distribution  or  superposition  of  sources  over  the  upper 
wing  surface.  Thus  the  velocity  potential  for  the  arrowhead  wing, 
irrespective  of  the  type  of  time  variation,  is 


0(x,y,+O,t)  =  - 


w(!  ,T|) 


an  as 


(8) 


where  w(|,t^)  represents  the  surface  variation  of  source  strength  and 


M(x  -  |)  PR 

h2  =  ap2 


R  =  ^x  -  |)^  -  P2(y  -  T])2 


The  region  of  integration  S,  in  equation  (8),  is  the  part  of  the  wing 
(in  It) -plane)  cut  out  by  the  reflected  Mach  cone  opening  upstream  of 
the  point  (x,y,+0)  as  shown  by  the  shaded  area  in  figure  2. 


For  the  present  problem  of  the  wing  performing  small  harmonic 
torsional  and  vertical  oscillations,  the  time  variation  of  source 
strength  is  given  by 

/ .  \  ia)t 
w(t)  =  e 

Thus,  as  shown  in  reference  1,  equation  (8)  takes  the  form 


$(x,y,+0,t) 


oM  ifiod 

2  2 
ap  Up 


where 
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and  the  values  of  w(l,q)  for  the  harmonic  case  from  equation  (5)  are: 
For  h, 

v(!,q)  =  icohQ  =  03  ho 

tr 

For  Ua, 


For 


«(= 


w(!,n)  =  UaQ 


W(5,t|)  =  — ^  (l  -  X^CCq 


For  certain  parts  of  the  region  of  the  wing  (shaded  in  fig.  2)  to 
which  equation  (9)  applies,  the  integration  can  he  carried  out  as  in 
the  case  of  the  infinite  wing  considered  in  reference  1.  However,  the 
entire  potential  0  does  not  appear  to  he  obtainable  in  terms  of  known 
functions.  Previous  mention  of  this  point  was  made  in  the  introduction. 
In  order  to  give  consistent  treatment  to  the  entire  region  of  integra¬ 
tion,  the  integrand  of  equation  (9)  may  he  expanded  in  Maclaurin's 
series  with  respect  to  cd,  with  the  result  that  each  term  of  the  expan¬ 
sion  can  he  readily  integrated.  Thus  equation  (9)  is  used  in  the 
following  form: 


0(x,y,+O,t) 


n=0 


(-if  (x  -  nm~2n 

(2n) !  (m  -  2n)  ! 


dTi  d| 


(10) 


where 


denotes  the  integral  part  of  m/2. 


Integrated  Form 

The  velocity  potential  to  the  third  power  of  cd  is  considered 
sufficient  for  a  large  number  of  practical  applications.  Thus  equa¬ 
tion  (lO)  to  the  third  power  of  .  cd,  with  any  cd  that  may  come 
from  w(  t  ,rj  ,t)  temporarily  neglected,  is 


0(x,y,+O,t) 


1 

R  + 


5 

al  R 


+  a 


2  T 


+ 


'3  R 


+  b0R  +  bf  R^dq  d£ 


(11) 
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ft 


where 


= 


_  ,  —  a)2  o  a?  3 

1  -  icox  -  x  +  i  x-3 


o  o>3  2 

=  ico  +  co^x  -  i  —  x 


a3  = 


.  *3 
-1  T 


ho  = 


_2  _3 

a>  .  or 
— ~  +  1  — ^  x 


2M‘ 


2M‘ 


\  = 


2M2 


The  question  may  he  raised  at  this  point  as  to  whether  equation  ('Ll) 
still  represents  a  solution  of  the  boundary- value  problem  under  discus¬ 
sion.  With  regard  to  this  question,  it  can  be  shown  that,  when  all  the 
terms  involving  cd  up  to  a  given  power  are  taken  into  account,  the  dif¬ 
ferential  equation  (l)  is  satisfied  to  this  power.  The  boundary  condi¬ 
tion  of  tangential  flow,  equation  (5),  is  satisfied  exactly,  regardless 
of  the  order  of  c5  considered. 


Putting  the  values  of  w(6,t))  given  in  equation  (9)  into  equa¬ 
tion  (ll)  and  performing  the  indicated  surface  integration  (over  the 
shaded  region  in  fig.  2)  yields  the  following  forms  for  the  velocity 
potentials  in  equation  (6): 


A)P 


„2  2 

-  3  y 


+  B-i  cos 


-1 


x  +  p2),y 
P(Xx  +  y) 


+  Br) 


-1 

COS 


x  -  pg\y 
p(\x  -  yT 


f  Fl(x,y) 


1  -  - 


Ua 


a 


fc  = 


n  Fi(x>y) 


> 


q2  2 

-  3  y 


+  D]_  cos 


-1  x  + 


P2*y 


P(Xx  +  y) 


+  D2  cos 


-1 


X  _ 

pTix 


-  32Ay 


71 


J 


(12) 
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The  coefficients  A,  B-j_,  B2,  C,  D-^  and  D2  are  given  in  appendix  A. 

In  order  to  perform  the  integrations  that  led  to  the  results  represented 
by  equation  (12),  the  use  of  characteristic  coordinates  (Mach  lines  in 
fig.  2)  was  found  convenient.  Reasons  for  using  these  coordinates  are 
given  in  appendix  A ,  where  a  brief  account  of  the  derivation  of  equa¬ 
tion  (12)  is  presented. 

In  the  limiting  case  where  the  Mach  lines  coincide  with  the  leading 
edges  of  the  wing,  that  is,  when  (31  =  1,  equations  (12)  take  the  fol¬ 
lowing  form: 

%=-l  A^x2  -  P2y2  I 


<P^  =  -  |(Cl  -  xoAi)\/x2  -  32y2 

J 

The  coefficients  A^  and  C-^  and  the  derivation  of  equation  (13)  are 
also  given  in  appendix  A. 


FORCES  AND  MOMENTS 


The  preceding  result  for  the  velocity  potential,  equation  (12) ,  is 
now  used  to  obtain  the  forces  and  moments  on  any  one  of  a  series  of 
wings  generally  referred  to  herein  as  arrowhead  type.  As  stated  in 
the  introduction,  the  wings  being  considered  are  those  that  may  be 
formed  from  the  delta  or  triangular  wing  by  cutting  the  rearward  part 
of  this  wing  in  such  a  manner  that  the  resulting  trailing  edges  lie 
ahead  of  the  Mach  cones  emanating  from  their  foremost  points.  Sketches 
of  several  wing  plan  forms  having  the  aforementioned  characteristics  are 
shown  in  figure  3* 

The  forces  and  moments  acting  on  a  wing  section,  such  as  section  y 
of  figure  1,  are  derived  in  general  form,  that  is,  in  the  form  applicable 
to  any  of  the  aforementioned  wings.  These  forces  and  moments  are  useful, 
for  example,  in  a  strip  flutter  analysis.  Total  forces  and  moments  are 
also  derived,  but  only  for  the  delta  wing.  The  damping  part  of  the 
total  moment  due  to  ocq  is  used  in  a  study  of  the  possible  loss  of 
aerodynamic  damping  (loss  indicates  that  single-degree  torsional  flutter 
is  possible)  of  an  oscillating  delta  wing,  which  infers  the  same  pos¬ 
sibility  with  regard  to  the  other  wings  considered. 
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In  deriving  the  expressions  for  the  section  and  total  force  and 
moment  coefficients,  it  is  convenient  to  employ  the  variables  p, 
and  pq  as  the  nondimens ional  quantities,  obtained  by  dividing  the 
variables  x,  y,  and  xq  by  the  chord  2b  of  the  midspan  section  of 

the  wing,  and  to  introduce  the  reduced  frequency  k  =  —jr. 


Section  Forces  and  Moments 


The  local  pressure  difference  between  the  upper  and  lower  surfaces 
of  the  wing  is  given  by  means  of  equation  (7)  as 

-  -MU +  h  t)  (l4 

The  expression  for  the  section  force,  positive  downward,  is  therefore 

m 


P  =  -2b 


Ap  dp 


and  the  section  moment,  positive  nose  up,  about  the  arbitrary  axis  of 
rotation  x  =  xq  is 

^  =  -4b2  P  1  Aj.  -  h0)Ap  d|!  (16) 

Jv/\  ' 

The  limit  of  integration  pp  in  equations  ( 15)  and  (l6)  has  the  fol¬ 
lowing  values  with  respect  to  the  different  plan  forms  shown  in  figure  3 


For  plan  form  A, 


For  plan  form  B, 


For  plan  form  C, 


For  plan  form  D, 


Pp  =  1 


^1  =  1  -  m 


-  V 

“1  ■  1  +  5 


ml  + 


l  v  _  ^2  (l  \  s 

2b  "  mg  f0r  m1  +  m^Sb  7  ^  v  ^  2b 
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After  equation  (12)  is  substituted  into  equations  ( 15)  and  (l6)  and 
the  indicated  integrations  are  performed,  the  results  can  he  written  as 


P  = 


.  2  2  io)t 

UpbU  k  e 


h0/T  \ 

/  \ 

—  (L1  +  1L^ 

+  a0(Lg  +  1L4J 

(18) 


and 


^  - 


4(»We“ 


■r(«l  +  iM2) 


(19) 


In  equations  (l8)  and  (19)  the  reduced  frequency  k,  on  which  the  Lpa 
and  Mj^'s  are  dependent,  is  related  to  o>  and  55  in  the  following 
manner : 


k 


to  _  bp^ 

u  V 


0} 


A  brief  account  of  the  derivation  of  and  expressions  for  the  force  and 
moment  coefficients  Lj  and  Mj_  (where  i  =  1,2, 3 A)  appearing  in 
equations  (l8)  and  (19)  is  given  in  appendix  B. 

In  equations  (l8)  and  (19)  the  coefficients  L-j_  4-  iL2  ,  +  iMg 

and  +  iL^  ,  +  iM^  are  the  section  lift  and  moment  coefficients 

associated  with  vertical  translation  and  rotation,  respectively.  The 
real  component  of  +  H2  ,  for  example,  is  in  phase  with  the  dis¬ 
placement  hg  and  the  imaginary  part  is  90°  out  of  phase  with  this 
displacement.  Similar  interpretations  can  be  given  to  the  remaining 
components.  The  imaginary  part  of  each  of  the  coefficients  is  propor¬ 
tional  to  the  aerodynamic  damping  force  or  moment  associated  with  the 
respective  motion. 


Total  Forces  and  Moments  for  Delta  Wing 

The  expression  for  the  total  force,  positive  downward,  on  plan 
form  A  of  figure  3  is 


P  =  -8b' 


v/X 


Ap  dp  dv 


0 


(20) 
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The  total  moment,  positive  nose  up,  about  the  arbitrary  axis  of  rota¬ 
tion  x  =  Xq  is 


M„  =  -l6b3 
a 


0  Jv/\ 


Ap  ^x  -  Xq)  dp  dv 


(21) 


After  equation  (12)  is  substituted  into  equations  (20)  and  (2l) 
and  the  indicated  integrations  are  performed,  the  results  can  be  -written  as 
(using  barred  letters  throughout  to  designate  total  forces  and  moments) 


P  =  - 8 pb^U^k2e 1CDt  ^(lx  +  iL2)  +  Oq(l3  +  iLjJ 


and 


% 


=  -8pb3U^k2eia>t|^^M1  +  iMg)  +  aQ^M3  +  iM^) 


(22) 


(23) 


The  derivation  of  the  total  force  and  moment  coefficients  Lj^  and  Mj_ 
(i  »  1,2,3, 4)  is  also  given  in  appendix  B. 


SOME  CALCULATIONS  AND  DISCUSSION 


It  may  be  of  interest  to  the  reader  to  examine  the  spanvise  dis¬ 
tributions  of  the  various  components  of  the  section  lift  and  moment 
coefficients  for  a  particular  case.  Thus,  the  components  of  equa¬ 
tions  (l8)  and  (l9)>  given  more  fully  in  equation  (Bl),  have  been 
evaluated  at  different  spanwise  positions  y  for  plan  form  A  of  fig¬ 
ure  3  for  the  following  set  of  conditions:  \  =  \fj~,  i±q  -  0-5, 

=  1.75,  and  k  =  0.04.  These  sample  results  are  plotted  as  functions 
of  the  span  position  in  figure  4.  The  spanwise  variations  of  the 
section- force  components  are  shown  in  figure  4(a)  and  the  corresponding 
variations  of  the  section-moment  coefficients  in  figure  4(b) . 

The  components  of  the  total  force  and  moment  coefficients  (given  by 
equation  (B4))  have  also  been  evaluated  for  the  previous  set  of  condi¬ 
tions.  The  results,  after  being  referred  to  the  nondimensional  span  of 
the  wing  2A,,  are  represented  by  dashed  lines  in  figure  4.  As  should 
be  the  case,  the  areas  under  the  dashed  lines  are  equal  to  those  under 
the  respective  section-component  curves. 

As  a  result  of  the  linear  dependence  of  the  total  force  and  moment 
coefficients  on  the  semiapex  angle  of  the  wing,  as  may  be  noted  in 
equation  (b4) ,  the  average  ordinates . of  the  distributions  given  in 
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figure  4  will  not  change  if  the  apex  angle  is  changed.  Only  the  shapes 
of  the  section  force  or  moment  distributions  will  be  affected  by  such  a 
variation  in  angle. 

Injthe  example  plotted,  the  total  component  of  moment  coeffi¬ 
cient  M 4  is  negative.  This  term  would  therefore  not  contribute  to 
the  aerodynamic  damping  but  would  act  as  a  source  of  energy  for  the 
oscillating  system.  This  result  is  significant  since  it  indicates  the 
possibility  of  single-degree  torsional  flutter. 

The  wing  plan  forms  discussed  in  this  paper  exhibit  the  possibility 
of  loss  of  aerodynamic  damping  in  torsion  for  certain  ranges  of  Mach  num¬ 
ber  M  and  axis-of-rotation  location  pq.  This  possibility  is  indi¬ 
cated,  as  mentioned  in  the  previous  paragraph,  by  the  negative  value  of 
the  torsional  damping-moment  component  M4  in  the  example  given.  Since 
a  wing  oscillates  as  a  rigid  body  at  a  fairly  low  frequency  of  oscilla¬ 
tion,  the  main  results  of  the  loss-of-damping  phenomenon  can  be  obtained 
by  maintaining  M4  to  the  order  of  l/k  only.  Thus,  there  is  obtained 
for  plan  form  A  of  figure  3  the  following  result  from  equation  (B4) : 


M4  = 


3P3k 


12(1^  -  l)u02  -  4|aQ(4M2  -  5)  +  3(2^  -  3)]  (24) 


_In  general,  torsional  stability  depends  on  the  sign  and  magnitude 
of  Positive  values  of  M4  indicate  stability  and  negative  values 

indi£ate  possible  instability.  The  borderline  case  is  thus  given 
by  M4  =  0. 

The  ranges  of  values  of  Mach  number  M  and  axis  of  rotation  Pq, 
for  which  M4  in  equation  (24)  vanishes,  are  shown  in  figure  5-  Only 
one  curve  is  obtained  for  all  delta  wings  with  supersonic  edges,  that 
is,  ^  1,  because  M4  is  dependent  on  X  for  magnitude  only,  as 
may  be  noted  in  equation  (24)  .  The  region  inside  the  curve  in  figure  5 
is  the  region  of  possible  instability.  For  convenience,  a  second 
ordinate  is  given  in  figure  5  showing  the  values  of  X  below  which  the 
leading  edges  become  subsonic.  With  the  aid  of  this  second  ordinate  it 
may  be  seen  that  a  delta  wing  with  a  semi  vertex  angle  of  45  is  the 
narrowest  triangular  wing  that  may  show  a  torsional  instability  and  yet 
retain  the  characteristic  (3X  >  1  because  the  uppermost  part  of  the 
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curve  in  the  figure  corresponds  to  a  value  of  X  =  1  when  (3 A.  =  1. 

For  a  narrower  wing  to  have  at  least  sonic  edges,  it  must  be  moving  at 
a  Mach  number  which  would  place  it  in  the  stable  region  of  figure  5 • 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  July  11,  1951 
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APPENDIX  A 

COEFFICIENTS  OF  VELOCITY-POTENTIAL  EQUATIONS 


A  derivation  of  equations  (12)  and  (13)  is  given  briefly  here. 
The  surface  integral  represented  by  equation  (ll)  when  written  in 
terms  of  the  characteristic  coordinates  u,v  becomes 


where 


u  =  Jjj(  I  -  Pn) 

uq  =  ^  (x  -  Py) 

v  =  H  (5  +  P^) 
vo  -  •§  (x  +  ^ 

The  limits  of  integration  in  equation  (Al),  from  the  characteristic 
coordinates  previously  defined  and  from  figure  2,  are: 

(A2) 
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where 

=  1  -  PX 

r  ~  1  +  px 

Putting  the  appropriate  values  of  w(u,v)  and  the  limits  of  integration 
given  in  equation  (A2)  into  equation  (Al)  and  performing  the  indicated 
integrations,  and  then  reverting  to  xy- coordinates,  yields  the  potentials 
in  the  forms  given  in  equation  (12).  Retaining  only  those  terms 
associated  with  h  or  a  that  involve  o>  up  to  the  third  power  yields 
the  following  values  for  the  coefficients  in  equation  (12): 


*  6MVL  '  J  6M2a2  ^ 

i  ^33p2(l8a2  +  3 Iff  +  15)  +  3cr(5a  +  3)]  - 

1  M%2f~p(38a  +  45)  +  3a(4a  +  9)] 


Bl(x,y)  =  -jTg 
cr 


(Xx  +  y)  -  i 


(XX  +  7)2  ‘  i^5/2  P(2CT  +  3)  +[  (A3) 


cl  (Xx  +  y)3  +  i  -XPpCd3/p[m2(2c  +  5)  +  3a  I  (Xx  +  y) 
J  48M2a>/^L  A 


B2(x,y)  =  B1(x,-y) 


1  +  3)  -  i 

3a  a^  24M^a:>  L 


Xx  .  Xx  55/ _  L  .  Xp  y  03  XxJa) 


a(a  +  3)  -  M2(26a  +  45)  +  9a 

-J  24m2  a^  — 


2  a  + 


19a  +.15) 


Continued  on  next  page 
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E>l(x,y)  = 


-  -jffrjiXx.  +  y)  -  ^y|(Xx  +  y)  +  --|y^(Xx  +  y)‘ 

i  Xx  +  y)2  +  i  ~—^(2a  +  3)(Xx  +  y)3  + 

«2..ts2  r  ~  ~~t 


-  ^-TA-[4M2(4g  +  5)  -  c(o  -  3)  (Xx  +  y)3  + 
48m‘ o'1  L  J 

+  5)  +  3oJ  (Xx  +  y)3  - 


^Sp(2  a  +  5)(Xx  +  y) 

48a  d 


Ko(x,y)  =  D,  (x,-y) 


where 


2  2 

a  =  p*X*  -  1 


For  the  case  in  which  the  velocity  components  normal  to  the 
leading  edges  are  sonic  or  in  which  the  Mach  lines  coincide  with  these 
edges  (that  is,  when  PX  =  1,  7  -  0,  or  cr  =  0) ,  it  may  be  noted  from 

equation  (A2)  and  figure  2  that  regions  II  and  III  no  longer  exist. 
Therefore,  the  potential  in  this  case  is  determined  by  integration 
over  region  I  only.  The  coefficients  of  equation  (13)  are  thus  found 
to  be: 


_2 

XO)  0) 


5  (x2  -  P2y2)  +  M2(tx2  +  2P2y2)J 


210M2 


7(x2  -  p2y2)  +  M^x2  +  2p2y2)J 


ci  -  If  ■ 1  i(4x2  -  p2y2>  -  ^[T(x2  - p2y2) +  ^  -  2p2y2J 


The  use  of  characteristic  coordinates  leads  to  very  simple  expressions 
for  the  limits  of  integration  given  in  equation  (A2) .  In  addition, 
their  use  simplifies  the  derivation  of  equation  (Ah) . 
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APPENDIX  B 


FORCE  AND  MOMENT  COEFFICIENTS 


Section  Force  and  Moment  Coefficients 


The  coefficients  Lj_  and  M-^  (where  i  =  1,2,3,  *0  are  obtained 
by  substituting  equation  (12)  into  equations  (15)  and  (l 6)  and  grouping 
the  results  in  the  form  given  in  equations  (l8)  and  ( 19) •  These 
coefficients  are  therefore  defined  as  follows: 


+  iLg  =  L]_ f  +  iLg 


L3  +  iLU  =  Ig'  +iL4*-  ±  +  2m0)(l1'.+  iV) 

Mj_  +  iMg  =  Mj_'  +  iM2*  -  2^^'  +  iL2*) 

M3  +  iM^  =  M3'  .+  iM^'  -  2|i0^L3*  +  iL^')  -  ^  +  2n0J  (Mj_  +  iMg) 

where 


Ll'  +  iL2'  =  ~h  2J^X  F1(2^2v)d(i  -  |  F1(2n1,2v) 


L3’  +  “V 


V  +  iM2’  = 


n  jj_i 

2/  F2(2u,2v)dn  -  i  F2(2|i1,2v) 

_  v/x  '  _ 

~  r^l  1 

4  J  F1(2n,2v)n  dn  -  i  2|a1F1(2(i1,2v) 


M3’  +  iM^1  =  -1  4  f1  F2(2h,2v)h  dji  -  |  2^(2^, 2V) 

LJvA  _ 

In  equations  (Bl)  the  quantity  Fp(2|i,2v)  is  obtained  from  the  quan¬ 
tity  Fp(x,y),  defined  in  connection  with  equation  (12),  by  replacing  x 
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by  2(1,  y 
cients  A, 


by  2V }  and  regarding  the  O)  associated  with  the  coeffi- 
and  B2  in  F-j_(x,y)  as  to*  where 


The  quantity  F2(2p,2V)  is  obtained  in  a  similar  manner  from  the 
expression 


F2(x,y)  =  C  ^x2  - 


the  right-hand  side  of  which  is  recognized  as  part  of  in 

equation  (12).  0 


Consider  now  the  following  general  definitions,  particular  forms 
of  which  are  contained  in  L1  and  (i  =  1,2, 3, 4): 


Emn  =  2“+n+l  v^n 


e  e+m+n 


F  =2 
rmn  * 


/  2  \e  m,  Nn  -1  Vi  +  P  a,v 

(e  xv)  03.  (Xm  +  v)  cos  p(Xlll  ;  v)  + 


n  -1  Up  +  P 


(-^)e  "l'h  -  »)“  -1 


G  =2 
mn 


m+n+2 


^  r  ^ 


v/x 


dp 


g  e  ^e+m+n+l 
mn  ~  ” 


A  M"(Xp  *  v)°  c°3"1  3,1  + 


(-32xv)=  Q  o"W  -  v)“  cos-1  do 


y  (B2) 


With  the  aid  of  equation  (B2),  the  1^'  and  '  (i  =  1,2,3, 4)  of 
equation  (Bl)  may  be  written  as 
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-2jrLx ' 

-2itL2  • 

-2jtL3  ’ 

-2  *V 

-2rtM1 • 

-2«M2’ 

-2nM3’ 

-2jtM4' 


~  (alE01  +  3lH0l°  +  P2F02°)  +  ^a2G02  +  a3G20  +  a4E03  + 
™"i  +  P3h03°  +  P4Fo4°)k2 


-PlF01  \  +  (alG01  "  a2E02  ■  a3E20  +  P2H02°  "  P3E03°)k 
71G01  +  72E02  +  73E20  +  B1H11  +  62H0]/  + 


0  0  0 
S3H02  +  54F12  +  85F03 


"(7; 

leoi 

4.  R  f  G  4.  ? 

rn  +  c 

>2F01i  +  83F02  )k  +  ( 

7 2G02  +  73G20  ‘ 

7,  E 

_ 

7  E  -1 

-  8.  H 

°  +  5  H  °  - 

■  8  ,F  ° 

-  8  F  1  - 

8„F 

no 

O 

5  21 

4  12  5  03 

6  13 

7  03 

8  o4 

ax| 

(E02 

"  G0l) 

0/0  0 
+  +  P2^f12  -  h02 

j]  +  [^2G03 

+ 

a3^2l 

a4^Eo4 

-  g03 

)  +  a5  ^E22  " 

G2l)  + 

0 

33H13  + 

P4(F14° 

-  H04°j 

k2 

Pl(H01 

-  F11°) 

K 

alG02  +  a2(G02  “  E03i 

)  +  a3(G20  - 

■  E2l) 

P2F12<~> 

'  P3(H03 

°.p 

13°  3k 

7igq2  + 

72^03 

“  G02 

)  +  73 (E21  - 

G2o)  + 

^1^21  +  SgH 

‘ii'!  + 

S3H12^ 

54(F22° 

-  h12°)  +  85 (f13° 

H  °’ 

H03  ; 

) 

7ll 

[G01 

“  E02) 

+  &1( 

H11°  -  F2i°) 

+  52(H011  -  FH^  + 

83  (e 

0 

[02 

°T 

-f12) 

I  + 

k 

—  \ 
72g03  +  73  G21  +  74^ 

/ 

^03  ~  eo4)  + 

r5(c 

'2.1  ' 

•  e22  + 

84®22  5^H]_3  + 

66(h13' 

D  „  0\ 

■  F23  J  + 

8t(b 

[03^ 

-  ^s1) 

+  &8(ho4°  -  Fl4°j 

k 
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where 


al  =  72 


[jftfc,  +  3)  +  a| 


«3  -  *  ') 


rft  a4 


- ^-g-3  [^^(18 a2  +  31a  +  15)  +  3a(5a  +  3)_ 


72 0  a 


72(3  V 


Pl  =?72 


P2  =  - 


45)  +  3a(4a  +  9|] 


IV" 

=- — £p5/a  f^2(2a  +  3)  +  a~l 


ft  M^X 

4  =  4i?7^ 


1^(20  +  5)  +  3a] 


7l  =  a 


r3  =  -  — 


M2).  ^  M2X 

72  - - 2~2  0  +  3)  74  =  -—4 

3P  a  24(3 


—  ^(2a2  +  190  +  15)  +  0(0  +  3)J 


75  =  ~24p 


|^3p(2  60  +  45)  +  9^ 
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81  =  '^72 


62  =  5X 


^  +  3) 


e  _  M2 
6  ' 


4m  (4<t  +  <§)  -  a(c 


-3^ 


sn.  f 

fvi  • 


c  p2x 


M£ 


s7  =  0 

'  48p2a 


^75  ^(So  *~i)  *  3jj 


68  “'S^75(20 

The  coefficients  in  equation  (B3),  that  is,  ot^,  (3-^,  7^,  and  5j_, 

originate,  respectively,  from  the  coefficients  A,  B-j_  or  B2,  C, 

and  D^_  or  D2  defined  in  appendix  A.  The  quantities  E^  and 

contained  in  equation  (B3)  are  easily  determined  from  equation  (B2) . 
The  quantities  and  H^,  however,  are  determined  by  integration, 

as  indicated  in  equation  (B2) •  The  integral  results  that  are  required 
in  equation  (B3)  are: 


64  =  - 


v£\ 

372 


2a 


G02  =  2 


3/2 


D2A,2  l  2 

+  p  V  - 


2^1  ^l2  -  P2v2) 

g03  *  i|(3^12  +  2P2v2)(^i2  -  P2v2'l/2 


cosh 


G20  =  8v2  klnfll-,2  -  P2V2)^  -  P2v2  cosh-1 


2k 
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TT  0  1  „  0  4a1' 2  2  ,  -1  U1 

H01  "  2X  F02  -  \  v  cosh 


n  0  __  -n  0  3gq_  __  .  .2v  ..  2  02..2 

h02  -  3x  f03 - 3  v  ri  "  p  v 


h03°  =  F04°  "  —  3m-x X2v2J|1;l2  -  P2v2  +  vV3<J  +  5)cogh_1  — ^ 


[o4&  =  5X  f05°  ‘  2,515A  _^(2a  +  5)  +  X3^ 


2v2  i/^i2  -  32v2 


t,  o  1  „  0  .  1  „  0  8a1/2  .,2  I  2  ^2^2 

H11  =  “^2  F03  +  2X  F12 - 3X~  V  FI  '  P  V 


tt  0  __  V  p  .0  l  .  -p  0  c-\2.J21]  .  [T|  2  p2A .2  , 

H12  “  _12X2  F04  +  3x  F13  "  ^2^L_X  V  ^  V  1  '  P  V 

v^(5o  +  3) cosh-1 


h13°  =-~2  f05°  +  II  Fl4°  ■ 


- 2  F05  +  la  Fl4 - 2 

20X  15X 


l6al/2 [T^  .  rtCw.k 


(22a  +  25) Xv4  + 


m3vV  iK2  -  p2v2 


0  10  10  10  4a1/2 

21  =  12X3  11 04  -  3X2  *13  +  2X  *22  -  3X2 


,2  2 

X  V  n1 


2  2 

p  v*  + 


V4(a  +  3)cosh" 


i  »v 

pv_ 


122  "sox3'05  6x2 
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-  --rtrlXv^  + 
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1  1 
01  2X 


FQ21  -  8 a1/2p2\v2  \fi x2  -  P2 


V 


H11  -  sc  F04  +  3  F13  “  2  f22 
2  1/2 


H, 


- (x2V2\xr  \p  -  P2v2  +  V^cosh”1 


1  i  F  i  l6pga1/'2v2 

03  =#x  F04 - 3 


i 


X2  +  Xv*(2cr  +Jl)  -  P2v2 


EKMTA 

A/p  l 


As  may  be  suspected  from  the  complicated  results  generically  represented 
by  Emn>  Fmn>  Gmn^  and  *tlie  section  lift  and  moment  coefficients 

are  not  simplified  to  any  large  extent  by  algebraic  combination.  The 
lift  and  moment  coefficients,  as  given  in  equation  (B3),  are  in  a  form 
which  is  perhaps  better  for  numerical  evaluation  than  a  combined  form. 


Some  cautionary  remarks  are  in  order  with  regard  to  equations  (Bl) 
and  (B3) .  In  evaluating  the  coefficients  and  Mj_  (where  i  =  3**0  > 

the  terms  involving  k^  should  be  ignored  because  the  lift  and  moment 
coefficients  associated  with  a  have  been  treated  completely  only  to 
the  same  order  as  those  associated  with  h/U.  In  equation  (B3)  only 
real  quantities  are  admitted.  Therefore,  the  following  restrictions 


must  be  noted  for 


v 

^1 


ii 


i 


2  „2  2  n 

m  -  p  v  =  o 


cos 


x  +  P  Xv 
3(Xmi  +  v) 


=  0 


Hx  -  P  XV 

COS  7T-rr - r-  =  It 

P(Xnx  -  v) 


1  h 

cosh  1  ^4=0 

Pv 
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Total  Force  and  Moment  Coefficients 


As  indicated  in  equations  (20)  and  (21),  the  total  lift  and  moment 
coefficients  Li  and  Mi  for  the  delta  wing  may  be  obtained  from  the 
previously  derived  Li  and  Mi  by  setting  ^  =  1  in  Li  and  Mi  and 

integrating  the  resulting  expressions  over  the  span  of  the  wing.  The 
results  of  the  spanvise  integration  combine  readily  to  yield  expressions 
which  are  much  simpler  than  the  corresponding  ones  in  the  section- 
coefficient  case.  Without  further  detail,  the  total  force  and  moment 
coefficients  are: 


L^  +  iL2  =  Li  +  iL2 


-\ 


L3  ♦  iL^  .  V  *  iV  -  (i  ♦  2„0)(Li  ♦  iL2') 

Mi  +  iM2  =  Mi  +  iMg '  -  2^^  '  +  1L2 

M3  +  mk  =  M3'  +  iMi/  -  2h0(l3’  +  iL^’J  "  ^  +  2no)(Mi  +  iM^ 


J 


where 


—  ! 

L3 


M 


1  = 


2X 

M2X  1 

333 

'  153?' 

X 

3k  ' 

235  k 

2 X 

333 

kx 

2M2X  1 

33k 

535 

X 

> 

33  " 

(B4) 
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k\  kl£\ 

3Pk  '  5p5 

1 6x 

1533 

2X  2M2\  , 

—  -  — —  k 


If  equations  (22)  and  (23)  had  been  put  in  the  conventional  form, 

Tj  =  C-^qS 

M  =  Cj^qSc 

where  S  is  the  area  of  the  wing  (for  the  triangle  S  =  4b  X)  ,  q  is 
the  dynamic  pressure,  and  c  is  the  wing  chord,  then  the  lift  and 
moment  coefficients  Cl  and  Cm,  obtained  from  equation  (B4) ,  would 
be  independent  of  X.  This  result  is  rather  remarkable  since  it  means 
that  the  lift  and  moment  coefficients  for  a  triangular  wing  with 
supersonic  edges  are  independent  of  the  vertex  angle  of  the  triangle. 
The  same  observation  was  previously  noted  in  reference  k. 
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